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Abstract 

We study chiral symmetry breaking using the standard gap equation, supplemented with the 
infrared-finite gluon propagator and ghost dressing function obtained from large-volume lattice 
simulations. One of the most important ingredients of this analysis is the non-abelian quark- 
gluon vertex, which controls the way the ghost sector enters into the gap equation. Specifically, 
this vertex introduces a numerically crucial dependence on the ghost dressing function and the 
quark-ghost scattering amplitude. This latter quantity satisfies its own, previously unexplored, 
dynamical equation, which may be decomposed into individual integral equations for its various 
form factors. In particular, the scalar form factor is obtained from an approximate version of the 
"one-loop dressed" integral equation, and its numerical impact turns out to be rather considerable. 
The detailed numerical analysis of the resulting gap equation reveals that the constituent quark 
mass obtained is about 300 MeV, while fermions in the adjoint representation acquire a mass in 
the range of (750-962) MeV. 

PACS numbers: 12.38.Lg, 12.38.Aw, 12.38.Gc 



* Electronic address: Arlene.Aguilar@ufabc.edu.br 



^Electronic address: | Joannis.Papavassiliou@uv.es 



I. INTRODUCTION 



The dynamical mechanism responsible for chiral symmetry breaking (CSB) in QCD has 
been the focal point of extensive research during several years Q-Q. The study of CSB 
in the continuum involves almost invariably some version of the Schwinger-Dyson equation 
(SDE) for the quark propagator (gap equation). This non-linear integral equation has a 
notoriously rich structure, being extremely sensitive to the details of its kernel; the latter 
is composed by the various non-perturbative ingredients entering into the gap equation, 
most notably the gluon propagator and the quark gluon vertex. As is well-known, the gap 
equation displays "critical" behavior: the support of the kernel throughout the entire range 
of integration must exceed a certain critical value in order to generate non-trivial solutions 
for the quark propagator 3] . Given that most of this support originates from the infrared 
region, i.e. around the QCD mass scale of a few hundred MeV, the study of CSB through the 
gap equation furnishes stringent probes on the various methods and models aiming towards 
a quantitative description of the non-perturbative sector of QCD. 

In recent years, a large number of independent large-volume lattice simulations have 
furnished highly non-trivial information on the infrared (IR) behavior of two fundamental 

Is theories, namely the (quenched) gluon and ghost propagators, 



ingredients of pure Yang-Mi 
for both SU(2) and SU(3) 
that (in the Landau gauge) the QCD gluon propagator and the ghost dressing function are 



16l-l2l|. In particular, these simulations have firmly established 



IR finite and non-vanishing 



22 



24]. 



Given that the lattice is expected to capture reliably the full non-perturbative information 
contained in the gluon and ghost propagators, it is natural to explore their consequences 
for CSB. To that end, in this article we use lattice results for these Green's functions as 
inputs for the gap equation, and study the emerging CSB pattern for quarks (fundamental 
representation) and for fermions in the adjoint representation. Specifically, we will employ 
the lattice data of 18], given that they perform SU(3) simulations for both the gluon and 
the ghost propagators. As it will become clear from the results presented in the main body 
of the paper, the analysis carried out here may be regarded as a serious test of the robustness 
of the aforementioned lattice results, and can serve as a characteristic example of the rich 
phenomenology that one may extract with them. 

The detailed implementation of the idea described above is far from straightforward, 
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mainly due to the complicated structure of the gap equation, which makes it difficult to 
determine its exact dependence on the aforementioned lattice ingredients (propagators) . In 
fact, of particular importance for the self-consistency of the whole picture is the role played 
by the ghost sector (see, e.g. ll|, and references therein). The way the ghost sector enters 
into the gap equation is through the fully-dressed quark-gluon vertex. Specifically, recall 
that, in virtually all treatments, the fully-dressed quark-gluon vertex is not obtained from 
the corresponding dynamical equation (the SDE of the vertex), but is rather expressed in 
terms of the quark propagator, such that its Slavnov- Taylor identity (STI) is automatically 
satisfied (this procedure is known as the "gauge technique" j^). The STI itself contains ex- 



plicit re 
kernel" 



erence to both the ghost dressing function and the so-called "quark-ghost scattering 



261 ]; the latter is given by its own dynamical equation, and, as we will see, its numer- 



ical impact to the solutions obtained from the gap equation is quite important. If, instead, 
one were to "abelianize" this part of the problem by assuming that the quark-gluon vertex 
satisfies a QED-like Ward identity rather than the correct (non-abelian) STI, the resulting 
gap equation would contain the gluon propagator as its sole ingredient, a fact that would 
lead to an apparent incompatibility, in the sense that the kernel would not exceed the critical 
value (no CSB), or would fail to generate realistic constituent quark masses jfj, 7, 10, 15|. 
The main results of the present article may be summarized as follows: 

(a) The non-abelian Ansatz for the full quark-gluon vertex depends explicitly on the 
the ghost-dressing function and the quark-ghost scattering kernel; the latter quantity has 
a rather complicated Dirac structure [271 ]. being composed by four independent form fac- 
tors [see Eq. (I2.10j) ]. As a consequence, the corresponding expressions for the form factors 
appearing in the Lorentz decomposition for the longitudinal part of the quark-gluon ver- 
tex [see Eq. (13. ip ] are modified (with respect to the case where the quark-ghost scattering 
kernel is set to its tree-level value); their full form is presented in Eq. (13. 5p . As a result, 
the gap equation acquires a more complicated structure, given in Eqs. ( 13.1ip -f l3~T2l) . To the 
best of our knowledge, both Eq. (I3.5P and Eqs. (I3.lip - fl3.12p appear for the first time in the 
literature. 

(b) The quark-ghost scattering kernel satisfies its own dynamical equation, which may 
be decomposed into individual integral equations for the various form factors entering into 
the gap equation. In turn, these integral equations depend, among other quantities, on the 
quark propagator, a fact which converts the full treatment of the problem into the solution 
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of a complicated system of various coupled integral equations. 

(c) In order to make the above system of equations more tractable, without compromising 
its main features, we retain only the dependence of the gap equation on the scalar form 
factor of the quark-ghost scattering kernel, discarding all other form factors. In addition, we 
choose a very particular kinematic configuration, which further simplifies the corresponding 
integral equation that determines the aforementioned quantity. The final "one-loop dressed" 
equation is given in Eq. A3. 16j) . and constitutes, to the best of our knowledge, a novel result. 
For the actual calculation of the scalar form factor we will use on the rhs of Eq. (13.161) the 



lattice results of 



18| , and then substitute the result (shown in Fig. [H]) into the gap equation. 



(d) A well-known endemic shortcoming of all approaches based on the gauge-technique 
is that the transverse (i.e. identically conserved) part of the (quark-gluon) vertex remains 
largely undetermined; this fact, in turn, distorts the cancellations of overlapping divergences, 
the multiplicative renormalizability of the Green's functions in question, and their compli- 
ance with the renormalization group (RG). The construction of the appropriate transverse 
piece has been carried out in detail for the case of QED 28], but no real progress has been 
made in a non-abelian context. As is common practice, the aforementioned problem is reme- 
died by multiplying (by hand) the kernel of the gap equation by an appropriate functions, 
which restores the desired properties. For the case at hand, the simplest quantity that ac- 
counts for the missing dynamics is the full ghost dressing-function. As we will explain in 
the corresponding subsection, this choice is dictated by the STI satisfied by the quark-gluon 
vertex, and enforces the correct RG behavior of the dynamical (running) mass obtained 
from the gap equation. It should be stressed that the inclusion of the dressing-function has 
a considerable numerical impact on the obtained CSB solution, boosting up the quark mass 
to phenomenologically acceptable values. 

(e) After substituting all necessary ingredients comprising its kernel, the resulting gap 
equation is finally solved numerically, for two different cases. First, we study fermions in 
the fundamental representation (quarks), obtaining a quark mass that in the IR is about 
300 MeV. Second, we consider CSB with fermions in the adjoint representation; the latter 



are particularly interesting, due to t 
and deconfinement they display [29 



re clear separation between chiral symmetry restoration 



31] . The corresponding mass obtained for the adjoint 



fermions is within the range (750-962) MeV, depending on the details of the quark-gluon 
vertex used. This values are not too far from what one would expect naively, given the 
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enhancement of 9/4 produced to the kernel of the adjoint gap equation due to the ratio of 
the Casimir eigenvalues of the two representations. 

The article is organized as follows. In Section [III we introduce the necessary notation, and 
review the general the structure of the gap equation. In Section II III we first construct an 
Ansatz for the quark-gluon vertex which makes full reference on the quark-ghost scattering 
kernel, and use this vertex to derive the corresponding gap equation. Then, the "one-loop 
dressed" approximation for the scalar part of the quark-ghost scattering kernel is set up, 
and the improvements necessary for restoring the correct RG properties are discussed in 
detail. In Section IIVI we present the main results of this work. In particular, after briefly 
reviewing the recent lattice results on the gluon and ghost propagators that enter into the 
gap equation, we proceed to the numerical solution of the gap equation. In Section [V] we 
discuss our results and comment on possible future directions. Finally, in an Appendix we 
analyze for completeness the structure of the gap equation within the framework of the pinch 
technique (PT) [32h34| or, equivalently, the background field method (BFM) [35 ]. 



II. GENERAL STRUCTURE OF THE GAP EQUATION 

In this section we will introduce the basic definitions and ingredients necessary for the 
study of the quark SDE (gap equation). Then, we will give a special emphasis in the con- 
struction of a general Ansatz for full fermion-gluon vertex, where its non-abelian character 
will be kept intact. Finally, we will incorporate it into the gap equation and explore its 
effects. 



A. Definitions and ingredients 

Let us first introduce the necessary notation. In covariant gauges, the inverse of the full 



quark propagator in the Minkowski space has the general form 



26] 



S- 1 {p)=j-m-E(p), (2.1) 

where m is the bare current quark mass, and S(p) the quark self-energy. It is common 
practice to decompose S(p), in terms of a Dirac vector component, A(p 2 ), and a scalar 
component, B(p 2 ), which allow us to define the dynamical quark mass function as being the 
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ratio Ai(p 2 ) = B(p 2 )/A(p 2 ), explicitly we have 



S-\p) = A^p 1 - Btf)I = A^ip 1 - M{p 2 n 



(2.2) 



where I is the identity matrix, and the term A~ 1 (p 2 ) is often referred to in the literature as 
the "fermion wave function". Note that, the fermion acquires a dynamical mass as long as 
B{jp 2 ) is different from zero. Therefore, the CSB will be signaled when we obtain B(p 2 ) ^ 0. 

a, n 




FIG. 1: The full fermion-gluon vertex. 

In addition, the gluon propagator A AtI/ (g) in the covariant renormalizable (R^) gauges, 
has the form 



P^(g)A(g 2 )+e 



q A 



where £ denotes the gauge-fixing parameter, and 



(2.3) 



(2.4) 



is the usual transverse projector. In this work, we are particularly interested in the Landau 
gauge which is reached when £ = 0. Moreover, the full ghost propagator D(q 2 ) and its 
dressing function F(q 2 ) are related by 



D(q 



, '/-"Or: 



(2.5) 



An essential ingredient in our study is the fermion-gluon vertex, represented in Fig. [TJ 
and given by 

T^Pi^Ps) = gT a T il (p l ,p 2 ,pz) , (2.6) 

with T a (a = 1, 2, ...,N 2 — 1) being the generators of the group SU(N) where the fermions 
are assigned. The matrices T a are hermitian and traceless, generating the closed algebra 



[T a ,T b ] = if abc T 



abcrpc 



(2.7) 
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FIG. 2: Diagrammatic representation of the fermion- ghost scattering kernel H(pi,p 2 ,p 3 ). 

where f abc are the (totally antisymmetric) structure constants. In the case of 577(3), and 
for fermions in the fundamental representations (quarks), we have that T a = X a /2, where 
A a are the Gell-Mann matrices. When fermions are in the adjoint, (T a ) bc = -if abc . 
The vertex T fJi (p 1 ,p 2 ,p 3 ) satisfies the fundamental STI 261 ] 

P3 r ^(Pi'P2,P3) = F(p 3 )[S- 1 (-p 1 )H(p 1 ,p 2l p 3 ) - H(p 2 ,p u p 3 )S~ 1 (p 2 )} , (2.8) 

where the fermion-ghost scattering kernel H(p\,p 2 ,p 3 ) is defined diagrammatically in Fig. |2} 
and is written as 

H a (p u p 2 ,p 3 ) = T a H( Pl ,p 2 ,p 3 ) . (2.9) 
The kernel H{m , p 2 , p 3 ) and the "conjugated" H(p 2 ,p 1 ,p 3 ) have the following Lorentz 



decomposition [27J (note the change p\ <-> p 2 in the arguments of the latter) 

H( Pl ,p 2 ,p 3 ) = X I + + X 2 j> 2 + X&^pl , 

H(P2,PuP 3 ) = X I - X 2 fa - X^ 2 + X 3 d^y 2 , (2.10) 

where the form factors Xi are functions of the momenta, Xi = Xi(pi,p 2 ,p 3 ), and we use the 
notation Xi(p,r,q) = Xi(r,p,q) and a^ u = ^[7^,7^] (Note the difference between a^ u and 
the usually defined a^ u = |[7 M , Jv])- 



B. The renormalized gap equation 

The SDE for the fermion propagator is diagrammatically represented in Fig. [3j Using 
the momenta flow and Lorentz indices indicated in Fig. [31 the gap equation can be written 

as 

S-\p)=p'-m-C r g 2 j r^S(k)T v (-p,k,q)A^(q), (2.11) 
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FIG. 3: The SDE for the fermion propagator given by Eq. (|2.1ip . The gray blobs represent the 
fully dressed gluon and quark propagators, while the black blob denotes the dressed fermion-gluon 
vertex. 

where q = p — k, j k = fi 2e (2Tr)~ d J d d k, with d = 4 — e the dimension of space-time, and rj? 
is the fermion-gluon vertex at tree level. C r is the Casimir eigenvalue of the given fermion 
representation (r = F for the fundamental, and r— A for the adjoint). More specifically, for 
the gauge group SU(3), we have Ca = 3 and Cp = 4/3. Note that, m is a current fermion 
mass, the same appearing in the QCD Lagrangian, and in case of m ^ in Eq. (12. lip , the 
chiral symmetry is explicitly broken. In this work, we will consider the case m = 0, i.e., the 
chiral symmetry is kept intact at the Langragian level. 

All quantities appearing in Eq. (12. lip are unrenormalized; they are related to their re- 



spective renormalized counterparts, denoted with a subscript "R", through the relations 36 1 

Sr(p;v) = Zp X (n)S(p), 

A R (q;ii) = Z^{^)A(q), 

F R (q;ii) = Z;\n)F{q), 

T R (p,k,q;fi) = Zi(fi)r 1/ (p, k, q) , 

g R (fi) = Z-\^g = Z^ZpZ^g, (2.12) 

where Zp, Za, Z c , Z\, and Z g are the corresponding renormalization constants. Substituting 
Eqs. (12~T2|) into Eq. OdT]) . we obtain 



fi* (P) 



Z F j> - Z l C I g 2 R [ ^S R (k)T Ru (-p, k, q)A%(q) . (2.13) 
Jk 

In addition, the STIs of Eq. (12.81) imposes the all-order constraint 

Z 1 = Z; 1 Z F Z H \ (2.14) 

where Zu is the renormalization constant needed for the quark-ghost kernel, i.e., 
H = Zjj l H R . Now, in the Landau gauge, both the quark self-energy and the quark-ghost 
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kernel are finite at one-loop; thus, at that order, Zp = Zh = 1, and, therefore, Z\ = Z^ 1 , 
i.e., at one- loop the quark-gluon vertex renormalizes as the inverse the ghost propagator. 
Imposing the above approximation in the Eq. (I2.13p . we obtain 



where we have suppressed the subscript "R" to avoid notation clutter. 

III. INFLUENCE OF THE GHOST SECTOR ON THE GAP EQUATION 

In this rather lengthy and technical section we study in detail how the ghost sector enters 
into the gap equation. To that end, in subsection IIII Al we will use the STI to determine 
the dependence of the form factors of the longitudinal part of the quark-gluon vertex on 
the corresponding form factors Xj appearing in the Dirac decomposition of the quark-ghost 
scattering kernel, given in Eq. (I2.10p . Then, we will use the resulting vertex in order to 
derive the most general expression for the gap equation, displaying the dependence on all 
form factors Xj. In subsection IIII Bl we derive the "one-loop dressed" expression for the 
scalar form factor X , which is the only one that will be considered in the ensuing analysis. 
The modifications introduced into some standard forms of the quark-gluon vertex vertex, 
and the form of the gap equations obtained with them are presented in subsection IIII CI 
Finally, the adjustments necessary in order to enforce the correct RG properties of the gap 
equation are discussed in subsection IIII Dt 

A. The full fermion-gluon vertex 

The most general Lorentz decomposition for the longitudinal part of the vertex 
r M (pi,p2,P3) can be written as |27j] 




r M (pi,P2,P3) = L i7 M + L 2 ($i -jfc)(pi -lb)? + L 3 (pi + L^a^ipx - p 2 f , (3.1) 

where Lj are the form factors, whose dependence on the momenta has been suppressed, 
in order to keep a compact notation i.e. Li = Li(pi,p 2 ,p 3 ). The tree level expression is 
recovered setting L\ = 1 and L 2 = L 3 = L 4 = 0; then, T^(pi,p 2 ,P3) = 7^. 



S-\p)=i>-Z- l C r g 2 [ ltl S{k)T u {-p,k 



(2.15) 



Jk 
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Due to the fact that the behavior of the vertex T fl (pi,p2,P3) is constrained by the STI 
of Eq. (12. 8p . the form factors Lj's appearing into the Eq. (13. ip will be given in terms of the 
form factors X;'s of Eq. f[2~TU|) . 

More specifically, using the standard decomposition of 5* _1 (p) expressed in Eq. (12. 2p . it 
is relatively straightforward to demonstrate that the rhs of Eq. (I2.8P becomes 

p£r> 1)P2 ,p3) = F(p 3 )[C I + Ojrfx + C 2 i> 2 + Oj^zM , (3.2) 

with 

C = -A( Pl ) (p 2 1 X 1 + PvP2 X 2 ) + A{p 2 ) (p 2 r Xi+p 1 -p 2 X 2 ) - B( Pl )X - B{p 2 )X ; 

d = A( Pl ) ( Pl -p 2 X 3 - X ) - p 2 2 A(p 2 )X 3 - B(p 1 )X l + B(p 2 )X 2 ; 

C 2 = A(p 2 ) ( Pl -p 2 X 3 - Xo) - plA( Pl )X 3 - B( Pl )X 2 - B{p 2 )X l ; 

C 3 = A(p 2 )X 2 - A( Pl )X 2 - B( Pl )X 3 + B(p 2 )X 3 . (3.3) 

On the other hand, contracting Eq. (13.11) with p 3 , we have 

p^ 3 T^p u p 2 ,p 3 ) = {pl-pl^l+M-pl^-L^-^ (3.4) 

Equating the right-hand sides of Eq. (I3.2p and Eq. (13 .4p . we can express the Lj's in terms 
of the functions A, B and X,'s. Specifically, 

Li = ^^^A(p 1 )[X + (pl- Pl - P2 )X 3 }+A{p 2 ){X + (pl-p 1 -p 2 )X s 

+ l^-^B(p 1 )(X 1 +X 2 ) + B(p 2 )(X 1 +X 2 ] 
L * = nr F } P3 \^ U(Pi)[(pI+Pi-P2)Xs - X ] - A(p 2 )[(p 2 2 + Pl -p 2 )X 3 -X ] 



2(p2-p?) 

+ 7 ^K,{B(p 1 )(X 2 -X 1 ) + B(p 2 )(X 1 -X 2 )y, 



m-pi 

L 3 = -£^L!a(p!) (plX x + Pl -p 2 X 2 ) - A(p 2 ) (p 2 2 Xx+px-p 2 X 2 ) + B( Pl )X - B(p 2 )Xo\ ; 
P2 Pi I ) 

U = ^-{a( Pi )X 2 - A(p 2 )X 2 + B( Pl )X 3 - B(p 2 )X 3 \ . (3.5) 



The standard approximation in the literature is to use the tree-level value of H(pi,p 2 ,p 3 ), 
which is equivalent to setting Xo = Xo = 1 and Xj = Xj = 0, for i < 1, in the above 
equation. In addition, the effects of the ghost-dressing are also neglected, by imposing 
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F(ps) = 1. In this limit, we obtain the following expressions 

_ A( Pl ) + A(p 2 ) 

u 



L 



3 



2 




A( P i) - 


^(p 2 ) 


2(P?- 


■pi) 


S(pi) - 


S(P2) 


Pi - 


pi 


0; 





L 4 = 0; (3.6) 



which leads to the so-called Ball-Chiu (BC) vertex 

r Bc(Pl'P2,P3) = g 7 

+ (Pl 2 ~ P f { KPi) - A(p 2 )\ + [S(p0 - £(p 2 )] ) . (3.7) 

Pi P2 I z J 

We will next insert into Eq. (I2.15P the general quark-gluon vertex of Eq. (13.1 j) with the 

expressions for the form factors Lj given in Eq. (I3.5p . Defining pi = — p, p 2 = k, and ps = q 

and taking appropriate traces, it is straightforward to derive the following expressions for 

the integral equations satisfied by A(p 2 ) and B(p 2 ) 

+ (3-8) 



Using that, for the gluon propagator in the Landau gauge, A£(q) = — z3A(g), and 
p M PyA M1/ (g) = —ih(p, k)A(q), where 

hM ^» 2 ~^\ ,3,0) 

we can cast the gap equation in the form 

/• 3p-(fc + p)L 4 -2(L 3 + L 4 )%,/c) . ) 

+ A ^(p)p-^(p) A(,?) ^ } j ' (3 - n) 
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B. The "one-loop dressed" approximation for Xq 

In order to evaluate Eqs. (13.1 ip and (13. 12[) further, it is necessary to determine the 
behavior of the form factors entering into the definition of the L, through the Eq. (13. 5p . 
From Eq. (I2.10p . one can projected out the form factors Xj in the following way 

Tr{H} 



X 

x 2 



4 ' 

plTrj^H} - ( Pl -p 2 )Tr{]/> 2 H} 

A\p\pl - (pi-p 2 ) 2 } 
plTi{pl 2 H} - (p.-p^Tijp'.H} 

A[p\pl - (pi-p 2 ) 2 ] 



Ti{a ap p°p p 2 H} 

3 mpi-PzY-PiPt)' 1 j 

It is clear from the diagrammatic representation given in Fig. [2] that H(pi,p 2 ,ps), and its 
form factors given in Eq. (13.131) . depend among other things, on the fully dressed quark 
propagator. Therefore, the treatment of the full gap equation given in Eqs. (13. lip and 
(I3.12p boils down to a complicated system of several coupled integral equations. In order to 
make the problem at hand technically more tractable, we will only retain one of the form 
factors given in Eq. (I3.13p . and study it in an approximate kinematic configuration, which 
simplifies the resulting structures considerably. Specifically, we will only consider the form 
factor X , and assume that pi = p 2 = p, and p = —p 3 /2, so that X = X . 

Note that the momentum p 3 coincides with the momentum of the gluon in the gap 
equation; therefore we will denote p 3 = q, and p = —q/2. To obtain a non-perturbative 
estimate for X , we will study the "one- loop dressed" contribution given by the diagram of 
Fig.H 

Denoting the full gluon-ghost vertex by Gff = 5 afc G^, the expression for X^ reads 

4' (P, P, q) = 1 " i (£) j ^ v {k)D{k - p) G v Tr{S(fc + q)T,} . (3.14) 
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FIG. 4: Diagrammatic representation of the quark-ghost scattering kernel, H(—q/2,—q/2,q), at 
one- loop. 

To evaluate this further, we will use the following approximations: (i) G v will be replaced 
by its tree-level value, G u = (k — p) u . Note that, since the full ^ u (k) is transverse (Landau 
gauge), only the p M part survives. (ii) For the vertex we will use a slightly modified 
Ansatz than that of Eq. (13. 7ft , in order to reduce the algebraic complexity of the 



resulting equation. 

Specifically, we will use as our starting point the Ansatz 

r„(g, -k-q,k)= 1 - ([A(k + q) + A(q)]^ + ^[A(k + q) - A(q)}(2^ + fl) , (3.15) 

which satisfies the STI of Eq. (12. 8p in the chiral limit, B = 0, and with the ghost sector 
switched off (F — H — 1). Since the second term on the rhs of Eq. (13. 15ft is proportional to 
k^, it vanishes when inserted into (13.141) . again due to the transversality of A /ii/ (/c), a fact 
that simplifies considerably the resulting expressions. 

According to the procedure discussed in the next section [see Eq. (13.221) ]. we will improve 
the Ansatz of Eq. (I3.15P by multiplying it by F(k), but keeping H = 1. Under these 
approximations, and after setting p = —q/2, Eq. (13. 14|) becomes (in Euclidean space) 



X U {q ) = l + £f[h-<<m A(k)D (* + I) F(k) 

J k L J 



A 2 (k + q)(k + q) 2 + B\k + q) ' 

(3.16) 



We will finally approximate A(k + q) — A(q) = 1, and B(k + q) — 0, then we obtain 

X^\q) = 1 + \c A g 2 q 2 J [I - f(k, q)}A(k)F(k)^^- , (3.17) 



where 



(3-18) 
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After carrying out the integral on the rhs of Eq. (13.17P one obtains an approximate expression 
for X^\q) in terms of A(fc) and F(k); the result will be reported in section HVl Note that if 
we had multiplied Eq. (13. 15 j) not only by F(q) but also by X^(q), as is done in Eq. (I3.22p . 
then, instead of simply computing the integral of Eq. (13. lTj) . one would have to deal with 
the more difficult task of solving an integral equation for the unknown X^\q). 



C. The gap equation with ghost-improved quark-gluon vertices 



There are two basic forms of the "abelianized" quark-gluon vertex usually employed in 
the literature: (i) the BC vertex, denoted by T^ c , whose closed form is given in Eq. ( 13. 7p . 
and (ii) the so-called Curtis and Pennington (CP) vertex |40|, to be denoted by T^ p . These 
two vertices differ by a transverse (automatically conserved term), namely 



rop(pi>P2,ps) = Ib (pi,p 2 ,p 3 ) + r£(p 1 ,p 2 ,p 3 ) 



(3.19) 



with 



Y(p 2 - pT) - (pi - pzYtyi +fa) 



[A{ P2 ) - A( Pl )} 



(3.20) 



where 



(p!-pI) 2 + 



pI + pI 



B 2 {p 2 ) | B\ Pl ) 
A 2 (p 2 ) A 2 ( Pl )_ 



(3.21) 



Evidently, under the approximations employed, the ghost effects due to F(p 3 ) and Xq\p s ) 
may be incorporated into these two vertices through simple multiplication of their form 
factors by F(p 3 )Xq\p 3 ). Denoting their "ghost-improved" versions by Tqy{puP2-,Pz) and 
^Qp{piiP2iPz)i respectively, we have 



rBc(Pi>P2,P3) = F(p 3 )X^\p 3 )T^ c {p 1 ,p 2 ,p- i ) , 
r£ P (pi,P2,P3) = F(p 3 )X [ Q ] (p 3 )T£ p (p ll p 2 ,p 3 ) . 



(3.22) 



Notice that, one recovers the vertex used in Ref. 
equations. 



ll| by setting Xq 1 '^) = 1 in the above 



Substituting the form factors of Pbc(P1)P2,P3) into Eqs. (13. lip and (13.121) . we arrive at 
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the following coupled system for A(p 2 ) and B(p 2 



A{p 2 ) = 1 + Crg'Z- 1 [ 

Jk 

B{p 2 ) = C r g 2 Z^ 



}C (p-k) 



A 2 (k 2 )k 2 + B 2 (k 2 \ 



}C B B c (k,p) 



(3.23) 
(3.24) 



A 2 (k 2 )k 2 + B 2 (k 2 ) 

where the kernel fCo(q) corresponds to the part that is not altered by the tensorial structure 
of the quark-gluon vertex, namely 

K {q) = A(q)F(q)X^ ] {q), (3.25) 



while the parts that are affected, JC BC (k,p) and /Cg C (fc,p), are given by [8| 



lC B A c (k,p) 



^[A(k 2 ) + A(p 2 )} [3p-k - 2h(p, k)} - 2B(k 2 )AB(k 2 ,p 2 )^± 



P 



-A(k 2 )AA(k 2 ,p 2 ' 



k 2 - 



(k-p) 2 k-p 



P 



+ 2-fh(p, k) 



P 



(3.26) 



)Cl c (k,p) = -B(k 2 )[A(k 2 ) + A(p 2 )] + 2 [B(k 2 )AA(k 2 ,p 2 ) - A(k 2 )AB(k 2 ,p 2 )] h(p, k) 



where h(p, k) is given by Eq. (13.101) and 

A(k 2 ) - A(p 2 ) 



AA(k 2 ,p 2 ' 



AB(k 2 ,p 2 " 



k 2 — p 2 



B{k 2 ) - B(p 2 
k 2 — p 2 



(3.27) 



Similarly, the effect of the vertex T^ P (pi } p2 } p3) is to replace the kernels K. BC (k,p) and 
/C| c (/c,p) appearing in Eq. f !3.23|) and (I3.24[) by /C^ p (/c,p) and JC B p (k,p), respectively, where 

K? A p {k,p) = ICr(k,p) + ^A(k 2 )AA(k 2 y^ 2 ^ !>J)1 



2p 2 



d(k,p) 



K c B p (k, P ) = JC B B c (k,p) + 3 -B(k 2 )AA(k 2 ,p 2 ) {k2 Pr 



d (k, P ) ' (3 ' 28) 
Note that the above equations are written in the Euclidean space. Specifically, the Wick 
rotation was performed by setting p 2 = —p 2 , A E (p^) = — A(— p\), A(— p%) = A E (— p 2 ), and 
B(—p 2 ) = B E (—p 2 ). Therefore, AA(k 2 ,p 2 ), AB(k 2 ,p 2 ), and h(p,k) change the sign under 
Wick rotation. For the integration measure, we used f k = i j kE , as a last step, we have 
suppressed the subscript "E" everywhere. 



D. Asymptotic behavior and renormalization group properties 

The study of the ultraviolet (UV) behavior of the dynamical quark mass, A4(p 2 ), pre- 
dicted by the coupled system formed by Eqs. (13.231) and (13.241) . reveals the need of one 
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final adjustment. Specifically, as is well known, the correct UV behavior for M.(p 2 ) is given 
by QQ 

M(x) 



c 

X 



111 ( A* 



(3.29) 



where c is a /i-independent constant, related to the chiral condensate (qq) u by [8] 



,8, 



47T 2 7j 



(qq). 



In 



A 2 



"7/ 



(3.30) 



and the mass anomalous dimension is given by 7/ = 9C T /(11C\ — 3/2C T rif). For the fun- 
damental representation, we have C r = 4/3 and then 7^ = 12/(11Ca — 2nj). Equivalently, 
we can rewrite it in terms of the first coefficient b = (HCa — 2rij)/(487r 2 ) of the QCD (3 
function, where 7/ = 3Cf/(167t 2 6). 

On the other hand, after setting Z\ — 1 and Xq 1 ' = 1, neglecting AA(k 2 ,p 2 ) and 
AB(k 2 ,p 2 ), and setting A(p 2 ) = A{k 2 ) — > 1, the asymptotic equation that the system ( 13"3 
f l3.24p yields for A4(p 2 ) is given, after the standard angular approximation, by 



M{p 2 



SC r g 2 

167T 2 



A(p 2 )F(p 2 ) / dk 2 M(k 2 )+ / dk 2 A{k 2 )F(k 2 )M{k 



(3.31) 



In order to we verify whether indeed Ai(p 2 ) of Eq. (I3.29|) satisfies Eq. (I3.3ip . we substitute 
M.{j) 2 ) into the rhs of Eq. (I3.3ip . together with the one-loop gluon propagator 

^2 



A- X (p 2 ^ 



p 



1 + 7i# In 



(3.32) 



where 71 
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3 Ca — C r nf)/327r 2 , and the one- loop ghost dressing function 



F-\p 2 ^ 



9C A g 2 



4 48tt 2 



In 



(3.33) 



After performing the above substitutions, it is straightforward to see that for asymptotic 
large values of p 2 , the dominant contribution comes from the first integral, and that the 
solution of Eq. f)3.3ip is indeed of the form given in (13 .29 p . but with one important difference: 
the anomalous dimension assumes the value 7/ = 36C r / (35Ca— 6C r n/), instead of the correct 
value given above. Note that the difference is due to the non-Abelian contributions (gluons); 
setting Ca = 0, the two values coincide. 

The reason for this discrepancy can be traced back to a typical ambiguity, intrinsic to 
the gauge-technique. Specifically, the standard procedure of constructing a vertex based 
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L leaves the transverse (i.e. auto- 



2a 



281 ] . While, in the presence of 



on the requirement that it should satisfy the correct S 
matically conserved) part of the vertex undetermined 
mass gaps, the infrared dynamics appear to be largely unaffected, this ambiguity is known 
to modify the ultraviolet properties of the SD equations [25]. Essentially, failing to provide 
the correct transverse part leads to the mishandling of overlapping divergences, which, in 
turn, compromises the multiplicative renormalizability of the resulting SD equations. The 



construction of the appropriate transverse part is technically complicated even for QED |40j . 
and its generalization to a non-abelian context is still pending. Given the expectation that 
the restoration of the transverse part should not affect the infrared dynamics, the usual 
short-cut employed in the literature is to account approximately for the missing pieces by 
modifying (by hand) the SDE in question. 

Specifically, for the case at hand, the correct asymptotic behavior for Ai(p 2 ) may be 
restored by carrying out in Eqs. (I3.23|) and (13.241) the replacement 

Z^Ka^p) lC A , B (k,p)F(p 2 ) . (3.34) 

given that, for large p 2 , the perturbative (one-loop) expression for F(p 2 ) is given by Eq. (l3.33|) . 

It is straightforward to verify that, with this modification, Eq. (I3.3ip yields the correct 
value for jf. Note that even though, strictly speaking, one needs to supply only the asymp- 
totic form given in (13.33!) . ^ * s natural to assume that the non-perturbative completion of 
this formula will be given by Eq. (14. 3p . namely the full F(p 2 ). 

At first sight it would seem that the modification introduced in Eq. f !3.34|) amounts to 
the counter-intuitive replacement Z~ x — > F(p 2 ), i.e. trading a (cutoff-dependent) constant 
for a (/i- dependent) function of p 2 [U[. Even though from the operational point of view this 
is indeed what happens, the idea behind is slightly more subtle. A more intuitive way to 
interpret Eq. (13.341) is that the corresponding kernels are modified due to the presence of the 
(unknown) transverse parts, whose additional contributions must be such that, when prop- 
erly renormalized, will effectively amount to the replacement given in Eq. (I3.34p . Needless 
to say, it would be very important to determine the precise mechanism that leads to these 
modifications, but this is at present beyond our powers. 

It is instructive to understand from a slightly different (but ultimately equivalent) point 
of view why the (minimal) factor that must be supplied to Eq. (I3.3ip is indeed F(p 2 ). The 
argument relies on the (expected) RG-invariance of the mass function A4(p 2 ). Specifically, 
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let us for the moment assume that the underlying theory is QED, and set F = 1 into 
Eq. (I3.3ip . Then, Ai(p 2 ) is RG-invariant, because, in QED, due to the abelian Ward 
identities, the product g 2 A appearing on the rhs of Eq. (13.311) is RG-invariant. In the case 
of QCD this is no longer true in general; instead, in the i?^-type of gauges the corresponding 
RG-invariant combination is given by g 2 AF 2 [ll|, . Therefore, the simplest way to convert 
the product g 2 AF into a RG-invariant combination is to multiply it by F, which restores 
the RG-invariance of Ai(p 2 ). 



IV. NUMERICAL ANALYSIS 



In this section we will first review some of the recent lattice data for the gluon prop- 
agator A(q) and the ghost dressing function F(q). Then, we will substitute them into 
Eq. ( 13.1 7p . in order to to obtain an estimate for xj^\q). With all necessary ingredients 
[A(g), F(q), andX^g)] available, we proceed to solve numerically the coupled system 
formed by the integral equations ( I3.23P and (I3.24p when the non-abelian modifications of 
the BC and CP vertices are implemented. We will solve the corresponding equations both 
for fermions in the fundamental (quarks) and in the adjoint representations; for the former 
case we will use the numerical results obtained to compute the pion decay constant and the 
quark condensate. 



A. Nonperturbative ingredients: gluon and ghost propagators from lattice and 
the form factor X^\q) 



In this subsection we will comment on the nonperturbative form of the three basic Green's 
functions entering into the gap equation, namely A(q), F(q), and xj^\q). 

We start by_ showing on the left panel of Fig. [5] the lattice data for gluon propagator 
obtained in 



18] . The lattice data presented there correspond to a SU (3) quenched lattice 
simulation, where A(q) is renormalized at fi — 4.3 GeV. In this plot, we clearly see the 
appearance of a plateau in the deep IR region. The IR finiteness of the gluon propagator has 



been long associated with the dynamical generation of an effective gluon mass 



22 



32 



42|. In 



fact, the above set of lattice data can be accurately fitted in terms of the following physically 
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FIG. 5: Lattice results for the gluon propagator, A(g), and ghost dressing, F(q), obtained in 
Ref. [l8| and renormalized at fi = 4.3 GeV. Left panel: The continuous line represents the gluon 
lattice data fitted by Eq. flU} when m = 520 MeV, g\ = 5.68, p x = 8.55 and, p 2 = 1.91. 
pane/: The lattice data for F(q) fitted by Eq. (g3j) using g\ = 8.57, m = 520 MeV, p 3 = 0.25 and, 
p 4 = 0.68. 



motivated expression [43 1 

A-V)=mV) + g 2 
with m 2 (q 2 ) given by 

2/ 2 

m (q 



UC A gl ^(q 2 + Pl m 2 (q 2 ) 



96vr 2 



li- 



ra 



q 2 + p 2 m 2 ' 



(4.1) 



(4.2) 



where the fitting parameters are m = 520 MeV, g\ = 5.68, p\ = 8.55 and, p 2 = 1.91. The 
parameter m acts as a physical mass scale, whose function is to regulate the perturbative 
RG logarithm; so, instead of diverging at the Landau pole, the logarithm saturates at a finite 



value 43|. In addition, for large values of q 2 , we recover the one-loop expression of the gluon 
propagator in the Landau gauge given by Eq. (13.321) . Note that, contrary to conventional 
masses, dynamically generated masses display a non-trivial dependence on the momentum 
transfer q 2 32] . In particular, m 2 (q 2 ) assumes a non-zero value in the III, and drops rapidly 



in the UV in a way consistent with the operator product expansion 

On the right panel of Fig. El we show the lattice results for F(q) obtained from [181 ]; 
evidently, F(q) saturates in the deep IR at the constant value [22|, [231]. The ghost dressing 



function is also renormalized at \i — 4.3 GeV, and the data can be accurately fitted by the 
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FIG. 6: Numerical result for the form factor Xq (q) of the quark-ghost scattering kernel given by 
Eq. (f3TT7D when a(p 2 ) = 0.295. 



expression 



F (q 



1 9 C K gl ln fq 2 + Pzm 2 (q 2 



(4.3) 



4 48vr 2 

where the dynamical mass is given by Eq. (I4.2p changing the parameter p 2 — > Pa. The fitting 

parameters are g\ = 8.57, m = 520 MeV, p$ = 0.25 and, P4 = 0.68. 

The last ingredient to be determined is the form factor xj^\q). We proceed substituting 

the fit for lattice data for A(g) and F(q) presented in Fig. [5] into Eq. (I3.17p . Then, for deter- 
mining the integral given by Eq. (I3.17p . we should fix the value of g 2 (p 2 ). We adopt the same 



procedure of Ref. 43], where it was found that the perturbative tail of effective coupling, 
determined from the lattice data, is compatible with four-loop perturbative calculation at 
MOM scheme, presented in j^. More specifically, we use a(p 2 ) = g 2 (p 2 )/An = 0.295. The 
numerical result for (q) is shown in the Fig. [HI and it can be fitted by 



a 2 + 



(4.4) 



1 + (q 2 /q 2 oT ' 

where oi = 1.0, a 2 = 1.27, q 2 = 0.027 GeV 2 and, p = 1.0 when q 2 < 0.202 GeV 2 ; while for 
q 2 > 0.202 GeV 2 , we choose a x = 1.40, a 2 = 1.0, g 2 = 0.45 GeV 2 and, p = 0.63. 

Xq 1 ^) shows a maximum located in the intermediate momentum region (around 
450 MeV), while in the UV and IR regions X^\q) — > 1. Although this peak is not very 
pronounced (notice the small scale in y-axis), we will see soon that it is essential for pro- 
viding to the kernel of the gap equation the enhancement required for the generation of 
phenomenologically compatible constituent quark masses. 
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FIG. 7: Left panel: The individual contribution of the ingredients composing ICo(q). The 
green dotted-dashed line represents the case where ICo(q) = A(g), the blue dotted line the case 
where ICo(q) = A(q)F(q), in the red dashed curve Ko(q) = A(q)F 2 (q) and, finally the black 
continuous line represents the case where /Co (q) assumes the full form used in our calculation 
lCo(q) = A(q)F 2 (q)x]j\q). Right panel: The corresponding dynamical quark mass generated when 
we use in Eqs. (|3.23p and f|3.24j) the different forms of ICo(q) are shown in the left panel. 

B. Chiral symmetry breaking in the fundamental representation and the pion 
decay constant 

We now proceed to the solution of the coupled system of integral equations denned by 
the Eqs.f l3.23j) and (I3.24p . Note that, in addition, and in accordance to the discussion given 
in subsection IIII D\ we carry out the substitution 



F(q) F 2 (q) 



(4.5) 



in the Eq. (13. 25ft . In other words, the kernel )Co(q) appearing on the rhs of Eqs. (13.251) and 
(I3.24p will assume the final form 

/C (g) = A(g)F 2 (g)X 1] (g). (4.6) 



Before solving the system formed by Eqs. (13.23]) and (I3.24p with IC (q) given by Eq. (I4.6p . 
it is instructive to study the numerical impact that each of the functions composing ICo(q) 
has on the value of the resulting quark mass. 

The results of this exercise are presented in Fig. [7J where on the left panel we show 
the support that JC (q) receives when we turn on one by one the Green's functions that 
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compose it. Without a doubt, the biggest numerical contribution comes from the ghost 
dressing function, F(q 2 ). Nonetheless, one should not underestimate the effect caused by 
the scattering kernel, xj^\q), which is responsible for a considerable contribution to the 
dynamical mass generation, as presented on the right panel of Fig. [7J On this panel, we 
show the corresponding dynamical quark masses that are obtained solving the system of 
Eqs.f l3.25j) and H3.24[) when fCo(q) assumes one of the four forms presented in the left panel. 
Notice that, neither JCo(q) = A(g) nor /Co(<?) = furnish the right amount of 

support necessary to trigger chiral symmetry breaking. The chiral symmetry is only broken 
when }Co(q) = A(q)F 2 (q) and, phenomeno logically compatible values are only obtained 
when the effects of xj^\q) are incorporated in the )Co(q). It is important to mention that, 
although xjj\q) does not provide a sizable support for JC (q) in the deep IR, the small 
contribution it furnishes in the intermediate region is enough for increasing the mass from 
Ai(0) = 190 MeV to Ai(0) = 295 MeV. This result is consistent with previous observations 
in the literature jsl, Q], stating that the support crucial for quark generation originates from 
the intermediate region of the integration momenta. 

Let us now return to the solution of the system formed by Eqs. (13.231 and (13. 24ft . 
Substituting A(g 2 ), F(q 2 ), and xjpfa) into Eqs.f l3T23|) and fl3T24"|) . with the modification 
of Eq. (I4.6p . and fixing the Casimir eigenvalue in the fundamental representation i.e. 
C r = Cp = 4/3, we determine numerically the unknown functions A(p 2 ) and B(p 2 ). 

On the left panel of Fig. [HI the red dashed line represents the numerical result for the 
quark wave function A~ 1 {p 2 ). As we can see, for large values of p 2 , the function A(p 2 ) goes 
to 1, in agreement with the discussion presented in the Section l"HI Dl In the opposite limit, 
A^ij) 2 ) develops a plateau saturating in a finite value around 0.85. 

The red dashed line in the right panel of Fig. [8] represents the corresponding dynamical 
quark mass Ai(p 2 ), obtained as the ratio B(p 2 )/A(p 2 ). One clearly sees that M.(p 2 ) freezes 
out and acquires a finite value in the IR, while in the UV it shows the expected perturbative 
behavior given by Eq. (I3.29f) and represented by the blue dotted-dashed curve. 

For the sake of comparison, we also solve numerically the coupled system when the non- 
abelian version of the CP vertex is employed. As we can see from the continuous black curve 
of Fig. [HI the qualitative behavior of the results obtained with the CP vertex is very similar 
to that obtained with the BC vertex. Note that both vertices generate phenomenologically 
acceptable values for A'l(O). More specifically, for the BC vertex we obtain A'l(O) = 294 MeV 
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FIG. 8: Left panel: The numerical solution for the quark wave function A~ 1 (p 2 ) in the fundamental 
representation when the non abelian versions of BC (red dashed curve) and CP (black continuous 
curve) vertices, are employed Right panel: The numerical solution for the dynamical quark mass 
A4(p 2 ). The red dashed curve represents the dynamical mass generated when the BC vertex is 
used, while the black continuous line is the solution found with the CP vertex. 



while for the CP vertex the value is slightly higher Ai(0) = 307 MeV. 

Once the behavior of the dynamical quark mass is determined, one may attempt to 
reproduce some of the phenomeno logical parameters that depend directly on it. Such a 
parameters is the pion decay constant f^, which measures the "strength" of the CSB. The 
pion decay constant is denned through the axial-vector transition amplitude for an on-shell 
pion, (0 |A° M (0)| Tc b (k)) — if ir k tl S ab . Making use of the method developed by Pagels and 



Stokar 



49[ , and Cornwall 50] , f„ can be expressed in terms of the dynamical quark mass as 



P 

J IX 



dy-. 



M{y) 



M{y) 



y dM{y) 
2 dy 



(4.7) 



An 2 J »[ y + M 2 (y)] 2 

Substituting in the above equation the numerical solutions for M.(p 2 ) presented in the 
Fig. |H1 we obtain f n = 64.3 MeV for the case of the BC vertex, while with the CP ver- 
tex we get /„- = 68 MeV. These values should be compared to the experimental value 



f n = 93 MeV [511 ] . Evidently, the obtained values underestimate f n by almost 30%. The ori- 
gin of this suppression could possibly be traced back to some of the approximation used for 
the quark-quark-pion proper vertex 50j. Some improved versions of Eq. (14. 7p can be found 



in the literature; in particular, we will employ the expression of 



52J , where a correction term 



23 



is added to Eq. (14. 7\i . More specifically, 

11 = # + 



(4.8) 



with 



5ft 



4vr 2 



dy 



dM{y) 
dy 



V 



l M\y) I 



(dM\y) 



dy 



y 2 M 2 (y) 



dM{y) 
dy 



(4.9) 



Adding this term to the expression of Eq. (14. 7p . we obtain /„. = 76.4 MeV for the case of 
the BC vertex, and /„- = 80.6 MeV for the CP vertex. Although the results are still below 
the experimental value, the correction added clearly contributes in the right direction. 
Finally, we will compute the quark condensate, which plays the role of the order parameter 

2 n 

for dynamical CSB. The quark condensate at scale of v = 1 GeV is given by |9[ 



(qq) (1 GeV 2 



dy 



yM(y) 



4tt 2 J A(y)[y + M : 



(4.10) 



Substituting again the solution M.{p 2 ) and A{p 2 ) presented in Fig. [8] into 
Eq. (14.101) . we obtain (qq) (1 GeV 2 ) = (211MeV) 3 when the BC vertex is employed and 
(qq) (1 GeV 2 ) = (217MeV) 3 . This value should be compared to the typical value of the 
quark condensate (qq) (1 GeV 2 ) = (229 ± 9MeV) 3 |53|. 



C. Chiral symmetry breaking in the adjoint 

Next we will solve the same system of integral equations formed by Eqs. f!3.23|) and (I3.24p 
for both BC and CP kernels, given by Eqs. (13.271) and (I3.28j) respectively, when the Casimir 
eigenvalue in the adjoint representation i.e. C T = C\ = 3. 

When one switches from the fundamental to the adjoint representation, the overall effect 
in the gap equation is an enhancement factor of 9/4 due to the difference in the correspond- 
ing Casimir eigenvalues. Of course, due to the nonlinear nature of the gap equation, the 
wave function A~^-(p 2 ) and the adjoint mass A^ a dj(p 2 ) are not obtained from their quark 
counterparts through simple multiplication by 9/4. 

The numerical results for the adjoint representation are shown in the Fig. [9J On the left 
panel, we compare the fermion wave functions, A^(p 2 ), when we use the modified BC and 
CP vertices. As expected, both solutions displays the right asymptotic behavior while in 
the IR limit the fermion wave functions saturate in smaller values compared to those of the 
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FIG. 9: Left panel: The numerical solution for the fermion wave function A^(p 2 ) in the adjoint 
representation when the modified BC vertex (dotted red curve) and the modified CP vertex (black 
continuous curve) are employed. Right panel: The numerical solution for the dynamical fermion 
mass M^p 2 ). The dotted red curve represents the dynamical mass generated when the BC vertex 
is used. The black continuous line is the solution obtained with the CP vertex. 

fundamental representation. In addition, we can notice that in the adjoint, A~^-(p 2 ) does 
not display anymore a minimum as it does in the fundamental representation. 

On the right panel, we show the fermion dynamical mass M&dj{p 2 )- We see that the 
infrared saturation of A^ a dj(p 2 ) occurs for higher values compared to the values of Ai(p 2 ) 
in the fundamental. In particular, when the modified BC vertex is employed, one obtains 
.M adj (0) = 750 MeV, while for the CP vertex A4 adj (0) = 962 MeV. Clearly we can see hat 
the latter values are higher than 9/AAi(0). In addition, it is interesting to notice that the 
adjoint representation is more sensitive to the change from the BC to the CP vertex, since 
the difference between the results obtained with the both vertices is much more pronounced 
here than in the fundamental representation. 



V. DISCUSSION AND CONCLUSIONS 



In this article we have studied the CSB pattern that emerges when the ga p eq uation is 
combined with the available lattice data for the gluon and ghost propagators 



Partic- 



ular attention has been payed to the way the ghost sector enters into the gap equation. In 
particular, a complete Ansatz for the longitudinal part of the quark-gluon vertex has been 
constructed, which, in addition to the ghost- dressing function [11], captures the full depen- 
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dence on the various form factors composing the quark-ghost scattering kernel. This new 
vertex has been used to derive a more complete version of the gap equation, containing ad- 
ditional contributions from the quark-ghost scattering kernel. This latter quantity satisfies 
its own dynamical equation, which is coupled to the gap equation in a complicated way. 

In order to reduce the complexity of the problem, we have restricted ourselves to the 
interplay of the gap equation with the only the scalar part, Xq, of the quark- ghost scattering 
kernel, deriving its "one-loop dressed" expression. The form factor Xq was then determined 
from this latter expression, for a special momentum configuration, using the lattice data for 
the gluon and ghost propagators appearing in it. The resulting expression for Xq, shown 
in Fig. El reaches its maximum at a momentum of a few hundred MeV, where it displays a 
25% enhancement compared to its tree-level value. We emphasize that the numerical effect 
of including Xq into the gap equation is rather sizable; indeed, as can be seen from Fig. [7J 
it accounts for about 30% of the final result for the quark mass. 

Finally, all ingredients are substituted into the gap equation, which is solved for the 
case of quarks (fundamental representation), and for fermions in the adjoint representation. 
The quark mass obtained is about 300 MeV, in good agreement with phenomenology, and 
is rather insensitive to the details of the quark-gluon vertex employed (BC vs CP). The 
corresponding mass obtained for the adjoint fermions displays a stronger dependence on the 
form of this vertex, varying between (750-962) MeV. 

It is important to emphasize that there is a crucial complementarity between using lat- 
tice data as input into the gap equation, and, at the same time, employing an Ansatz for 
the vertex that captures suitably the dependence on quantities such as the ghost-dressing 
function and the quark-ghost scattering kernel. Indeed, the additional dependence on the 
ghost sector stemming from the vertex would be insufficient for getting the correct CSB if 
one were not to use the lattice results, which capture fully the underlying dynamics, dis- 
playing a sizable enhancement with respect to other non-perturbative approaches. Similarly, 
no realistic CSB pattern can be obtained if one were to substitute the lattice ingredients 
into the gap equation obtained from a less elaborate quark-gluon vertex, i.e., one that fails 
to include the effects due to the quark-ghost scattering kernel. Thus, at least within our 
framework, it is the interplay between these two points that finally provides an adequate 
description of CSB. 

Amplifying on the previous point, let us mention that the SDE analysis presented in 22], 
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even though it reproduces qualitatively the lattice results, and most importantly accounts 
for the observed IR finiteness of the gluon propagator and the ghost- dressing function, 
it underestimates the size of both quantities by a significant amount. Clearly, it is an 



important challenge for the SDE approach of the PT-BFM formalism [22|, |42|, [54J to eliminate 
the aforementioned quantitative discrepancy. Perhaps the most obvious step would be to 
extend the analysis of 22j beyond the "one-loop dressed" approximation, given that there 
is no a-priori guarantee that the omitted "two-loop dressed" contributions are numerically 
depreciable. Needless to say, from the technical point of view, such an attempt would 
constitute a formidable task. 

Another possible source of enhancement for the Green's functions in question may be 
related to the non-trivial structure of the vacuum, and in particular with the presence of 



es. These classic field configurations are 
and are known to affect 



56| and CSB 



solitonic structures, such as vortices or monopo 
closely linked to the mechanisms of confinement 
the shape and size of the fundamental Green's functions of the theory 58]. A particularly 
instructive example of how to include such effects at the level of the SDEs has been presented 
in [59]. 



Appendix A: The gap equation in the PT-BFM language 

Given the recent reformulation of the SD series within the PT-BFM framework 
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601 ] . it is conceptually interesting to re-express the gap equation using the PT-BFM 
terminology. Let us stress from the beginning, however, that, unlike what happens in the 
case of the gluon self-energy 54j, where the corresponding SD equation in the PT-BFM 
formalism is vastly different from the conventional one, the gap equations obtained within 
both formalisms are completely identical. 

Let us start the discussion by pointing out that the PT quark-self energy coincides with 
the conventional one in the Feynman gauge, both perturbatively (to all orders), as well as 
nonperturbatively. The reason has been explained in detail in the related literature; here it 
should suffice to mention that in the Feynman gauge there are no pinching momenta, and all 
three-gluon vertices appearing in the quark-self energy are "internal" , in the sense that all 
legs are irrigated by virtual momenta, and, therefore, they should not undergo the standard 
PT decomposition, a key ingredient in the construction of the PT gluon self-energy. An 
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equivalent way of saying this in the BFM language 35] is that, unlike gauge fields, fermionic 
fields are not split into a background and a quantum component. Therefore, the BFM 
fermion propagators are the same as the conventional ones (in all gauges). However, given 
that all nonperturbative ingredients we will use come from lattice simulations in the Landau 
gauge, the gap equation we study here is not the genuine PT gap equation. Away from the 
Feynman gauge, one must switch to the BFM language, or, equivalently, to the generalized 
PT. In any case, the central result remains the same: the gap equation in the Landau gauge 
is the same as the conventional one. 

The PT-BFM gluon self-energy, denoted by A, behaves in many aspects as that of QED; 
in particular, the product g 2 A is RG-invariant (for any choice of the BFM gauge-fixing 
parameter). Of course, the propagator appearing in Eq. ( 12.15)) is not A but rather the 
conventional A; indeed, background field propagators do not propagate inside loops, only 
quantum ones. However, there exists a set of powerful identities that allows one to establish 



some important connections 41 



61 



621 ] . Specifically, 

A(g 2 ) = [l + G(g 2 )] 2 A(g 2 ) ; 
F-\q 2 ) = 1 + G{q 2 ) + L{q 2 ) . 



(Al) 



The functions G(q 2 ) and L(q 2 ) are the two form factors of a particular two-point function, 



denoted by A^ u (q), defined as [54j . 
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KM) 



-iff 



''C A [ H^D(k + q)A^(k)H m/ (k,q) 
Jk 



9,MQ 2 ) + qj ^L(q 2 ^ 



(A2) 



Turns out that the function L(q 2 ) is sub leading both in the IR and in the UV; therefore, 
for the purposes of this argument can be safely neglected 41] . Then, the combination of the 
two identities in (lAlj) leads to the approximate relation 



A(q 2 ) = A(q 2 )F 2 (q 2 ^ 



(A3) 



leading to exactly the same conclusion as before. 

Finally, the above discussion may be recast in the more intuitive language of an effec- 
tive (running charge), traditionally employed in QED-inspired studies of QCD. From the 
(dimensionful) RG-invariant quantity g 2 A introduced above, one may define a nonpertur- 
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bative effective charge, denoted by a(g 2 ), as 

2 a / 2n _ 47T«(g 2 ) 



(A4) 



q 2 + m 2 (q 2 ) ' 

where m 2 (q 2 ) is the momentum-dependent dynamical gluon mass, whose value in the deep 
IR is about (500 - 700) MeV. Using Eq. flMJ), and after implementing Eq. (pTMj) . the 
alternative (and completely equivalent) form of Eq. (13.231) reads p, [7| (setting q = p — k) 

a(q 2 ) fC A (k,p) 



A(p z 



1 + AnC r 



(AS) 



If. q 2 + m 2 (q 2 ) A 2 (k 2 )k 2 + B 2 (k 2 ) ' 

with an exactly analogous expression for B(p 2 ). Note that finally there is no explicit ref- 
erence to F(q 2 ), because it has all been absorbed into the definition of the effective charge 
a(q 2 ). The RG-invariance of this equation can be easily established, given that both a(q 2 ), 
the gluon mass m 2 (q 2 ), are RG-invariant. The final inclusion of the Xq^ into Eq. (1A5j) is 
straightforward. 
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